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In this letter we study the effect of time-reversal symmetric impurities on the Josephson super¬ 
current through two dimensional helical metals such as on topological insulator surface state. We 
show that contrary to the usual superconducting-normal metal-superconducting junctions, the sup¬ 
pression of supercurrent in superconducting-helical metal-superconducting junction is mainly due 
to fluctuations of impurities in the junctions. Our results, which is a condensed matter realization 
of a part of the MSW effect for neutrinos, shows that the relationship between normal state conduc¬ 
tance and critical current of Josephson junctions is significantly modified for Josephson junctions 
on the surface of topological insulators. We also study the temperature-dependence of supercurrent 
and present a two fluid model which can explain some of recent experimental results in Josephson 
junctions on the edge of topological insulators. 


Critical currents in the mesoscopic superconducting- 
normal-superconducting (SNS) Josephson junctions have 
been studied widely [l| as an important property which 
can identify important physical properties of SNS Joseph¬ 
son junctions. Different regimes of coherent length, mean 
free path in the normal region and length of the normal 
region have been considered both experimentally and the¬ 
oretically. One of the important results in this regard is 
the relationship between the normal state conductance 
and the critical current of the Josephson junctions @11 
which is well established both theoretically and experi¬ 
mentally. The discovery of topological insulators (TIs) 
and the novel electronic band structure on their surface 
have led to many investigations on the unique features of 
electronic transport, such as quantum anti-localization, 
in this new class of materials @-@|. The experimental 
study of transport properties of the TI edge states has 
been quite challenging though as the bulk carriers con¬ 
tribute significantly to the transport in most of the exper¬ 
imentally accessible topological insulators. Superconduc¬ 
tivity in TIs have been another widely studied area of re¬ 
search. The realization of superconductivity in doped TIs 
@[I3, as well as the possibility of making heterostruc¬ 
tures of TIs with superconductors Els are important 
developments in this regard. Theoretical predictions such 
as topological superconductivity in doped TIs (l^ and 
presence of Majorana zero-mod ellJl in Josephson junc¬ 
tions through TI surface states [ISl . Il6| motivated many 
experimental studies. 

The presence of helical edge states on the boundary of 
TIs is by now well established. The surface state band 
structure which resembles massless fermions, provides a 
unique platform to realize phenomena previously stud¬ 
ied in high energy physics, such as axionsfl^ and super- 
symmetry [3 [l9| . in much more easily realizable con¬ 
densed matter system. An important property of such 
states is the absence of back-scattering from electrostatic 
impurities, which is enforced by the strong correlation 
between spin and momentum. A reversal of momentum 
needs a reversal of spin and since eleetrostatic fields can¬ 


not flip the spin, this helicity conservation forbids back- 
scattering. This, in turn, makes electronic transport 
through surface states of TIs insensitive to non-magnetic 
impurities. Given this situation, it might seem that such 
impurities will not affect the supercurrent through sur¬ 
face states of TIs as well. In this letter we show that con¬ 
trary to the situation for the normal state, non-magnetic 
impurities will affect the supercurrent carried by the sur¬ 
face states of TIs. This dynamical effect resembles yet an¬ 
other phenomenon familiar from high-energy physics in 
the context of neutrino oscillations known as Mikheyev- 
Smirnov-Wolfenstein (MSW) effect in topological insula¬ 
tors. We will show that the fluctuations of the impurities 
lead to a a renormalization of the Fermi velocity. This in 
turn means that the optical length of the junction (de¬ 
fined via the phase of the wave function) is larger than 
the geometric length. The modification of the phase also 
modifies (via matching conditions) the energy eigenval¬ 
ues of the Andreev states. This is the essence of our 
result. In fact it has been noticed before that oscilla¬ 
tions of impurities might affect the critical temperature 
in superconductors!^, but effect studied here as the sole 
mechanism for impurities to change the supercurrent in 
helical metals has not been considered before. Our results 
can be used to interpret the measurements on TI Joseph¬ 
son junctions which are currently the focus of many ex¬ 
perimental studies (2 iI - E^ . 

The low-energy effective Hamiltonian of TI surface 
states reads as 

TLs = vf cr ■ k (1) 

where cr = are the Pauli matrices in the basis 

with ipa being the electronic state with spin a 
localized on the surface of the TI. The low-energy effec¬ 
tive Hamiltonian describing a Josephson junction on the 
surface of the TI, with supercurrent along x, is given by 

H = {-ivpV ■ a - h)t3 + An{x)Ti + Aj{x)t2, (2) 
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taken to act on the electronic state of the form 
ipi, 'i/'j, ■ Here Afiix) and A/(a;) are real and 

imaginary parts, respectively, of the induced supercon¬ 
ducting gap A = Ar + iAj. In ([ 1 ]) and Q, vp is the 
Fermi velocity. As for the matrix structure, ai act on 
physical spin space whereas the act on the supercon¬ 
ducting particle-hole space. As the Hamiltonian in ([2]) 
is invariant under translation along y, the momentum 
ky in this direction is conserved. The low-energy An¬ 
dreev states in the junction thus correspond to ky = 0 
and kx close to the two fermi wave vectors, kx = for 
CTa, = 1 and kx = for = —1. Notice that since 
Gx commutes with Hamiltonian ([5]) we can decouple the 
low-energy effective Hamiltonians into two independent 
sectors corresponding to electron-hole states close to right 
or left fermi points. Here we will focus on one of the ef¬ 
fective Hamiltonians, but the other independent one can 
be similarly studied. Since we are aiming for the effect 
of temporal fluctuations it is more efficient to use the 
corresponding 1 -f- 1 -dimensional action given by 


—V(x — oo) $ (x, t) ir 2 <h (x, t) 

Sint = J dxdt dxV{x - ao) $ (x,t) 1 x 2 $ (x,t) 

Sosc = \mi j dt ^(t) e(t) (4) 

Since we are interested in localized impurities, both the 
impurity potential V{x — a) and the resulting electric 
field E{x) = dxV{x — a) are localized in space and the 
dynamics of the impurities are captured by the harmonic 
oscillator action Sosc (AI/ and ui are the mass and the 
harmonic oscillation frequency of the impurities). Cou¬ 
pling of Andreev states and the harmonic oscillations of 
impurities are captured in Smt and will lead to a self¬ 
energy correction in So (see supplementary materials), 
leading to an effective action of the form 


S= dxdt 


(1 -I- Ei(x)) $ (iTidt + M (x)^ 


$ 


+ $X29x$ + $^(x)S 2 (x) $ (5) 


5 ' = 




ITiDt + T2Dx + M {x) $ (x, t) 


(3) 


where $^(x,t) = {x,t ),and (j)a{x,t) = 


ipa {x,t) with a =t, 4.5 are the fermionic field op¬ 

erators for excitations close to the right or left Fermi 
points. Dfx = dy — ie is the covariant derivative. 
Notice that ([31) is of the standard form of a Dirac action 
+ M)$ if we identify 7 ° = xi, 7 ^ = — 1 x 2 , 7 ® = 
X 3 . Further, in ([S]), $(x,t) = $'l'(x,t) 7 ° = $'l'(x,t)xi 
and M (x) = A/{(x)xo -I- iA/(x)x 3 . The Fermi velocity 
has been set to 1 by scaling x, or equivalently, the mo¬ 
mentum kx- The effect of charged impurities is captured 
by Ao = V{x — a(t)) where a(t) identifies instantaneous 
position of the impurity. As we will see below, in order 
to capture the effect of impurities on the supercurrent 
in the junction, we should consider the natural fluctu¬ 
ations in the position of the impurity. For small fluc¬ 
tuations a{t) = ug -l- ^(t), the impurity potential reads 
as V{x — a{t)) « V{x — Oq) -I- dxV{x — ao)i{t). As 
we will show below, the impurities can only affect the 
supercurrent as a result of their temporal fluctuations. 
We would like to note that such treatment of impuri¬ 
ties and its effect on superconductivity in normal metals 
have been considered long before 113. But here we show 
that as a result of helical band structure of the surface 
states of the TIs, the temporal fluctuations are the sole 
mechanism through which impurities can affect the su¬ 
percurrent. The action, including the effect of fluctuating 
impurities, read as 


S — Sq ~\- Sint + So 


So = 


J dxdt $ (x, t) (j,Tidt + T 2 dx + M (x)^ $ (x, t) 


dxdt V (x 


ao) $ ix2 $ 


To the lowest nontrivial order in perturbation theory, the 
self-energies can be calculated as Ei(x) ~ ^ 


S2(x) 


e^E^{x) 

2'KIjJ'^Mi 


R 


log(t^) - 1 


where R is the length 


scale over which E{x) = —dxV{x — a) is non-zero (see 
supplementary materials). What is important for us is 
not so much the specific formulae for these self-energies, 
but that the general form of the effective action is as given 
in (©, with the self-energies as corrections concentrated 
around the impurities. 

Notice that the modification due to Ei does not af¬ 
fect the spatial derivative term for the electrons. This is 
because the oscillator variable ^(t) does not have a spa¬ 
tial dependence. Thus although the free fermion action 
$(* 7 ^ 9 ^-|-M)$ has a Lorentz-type symmetry (albeit with 
vp in place of the speed of light c) the interactions with 
impurities, and hence corrections, do not respect this 
symmetry. The temporal and spatial derivative terms 
can be renormalized differently, and then with a possible 
scaling of the field, the action can be brought to the form 
The general result is that the effect of fluctuating im¬ 
purities appear as renormalization of the effective Fermi 
velocity and of the size of superconducting gap in the 
region where the electric field of the impurity potentials 
are present. The change of Fermi velocity can be viewed 
as a “refractive index" for the electron. There is then an 
additional phase change acquired for the wave functions 
as the effective “optical length" is modified. This will be 
the essence of how the Andreev states are modified. 

It is well known that for massless particles, propagat¬ 
ing at the speed of light in vacuum, the primary effect of 
interactions is to generate a refractive index rather than a 
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mass (which is usually forbidden for symmetry reasons). 
A similar situation is obtained even for massive particles 
in the ultrarelativistic limit. Our argument is that, for 
the surface states in a TI which have a Lorentzian sym¬ 
metry (with c —>■ vp), again a refractive index is precisely 
what we should expect as the primary effect of interac¬ 
tions. 

The situation here is closely related, conceptually, to 
how interactions with matter modify neutrino oscillations 
as in the MSW effect There are two ingredients 

to this. First, the neutrinos acquire a refractive index 
which can be calculated in terms of the forward scattering 
amplitude, a calculation which mirrors, mutatis mutan¬ 
dis, what is given in the supplementary material. The 
refractive index for electron-neutrinos is different from 
that for other flavors due to charged current interactions 
with matter. The resulting difference in the phase of the 
wave functions modifies the neutrino mass-eigenstates in 
matter, and hence the oscillations between different fla¬ 
vors of neutrinos. The second part is a resonance effect 
which can enhance the mixing of flavors in matter, even 
up to the maximal mixing. For us, we have only one fla¬ 
vor to consider, so the situation is simpler; there is no 
resonance part, but the phase enhancement due to the 
refractive index is similar and can modify the matching 
conditions (and energies) for the Andreev states. 

We will now use the effective action ([S]) to study the 
effect of impurities on the supercurrent. The Andreev 
eigenstates are determined by the effective Schrodinger 
equation which follows from (O, 


The eigenstates can be expressed in each region as the 
superposition of ^'^(a:) as 'f'i(a:) = 

where f = I, II and III, corresponding to the three regions. 
We define the transfer matrices: 


fBiE,A) = 




\/E -\- Kv —A*/^E -I- Kv 
—Ajy/E Kv \/E KV 


fnix) = 


0 


0 


( 11 ) 


fs{x) 


,IKX Q 

0 e-*' 


where k = y/ E'^ — Aq jv is the wave vector in the su¬ 
perconducting region, = E jg the averaged 

wave vector in the normal (TI) region of the junction 
and 4>i{x) = is the phase resulting from 

the static impurity. In terms of these matrices, the eigen¬ 
states in the three regions can be written as 


Mx)=fB{E,Ao)fs{x) 


■ Ai' 


d'll(x) = fn{x) 


’A"' 


( 12 ) 

(13) 


^m{x)=fB{E,AoE'^°)tix-Xp) 


- ^iii ■ 
^iii 


(14) 


The boundary conditions which determine the spectrum 
of the states in the junction are 4/j(0“) = 4'n(0+) and 
4'n(a;“) = 4'in(a;p). These continuity equations reduce 
to 


-[{ivdx + V{x - ao)) T 3 -h Apn -h A/T 2 ]'i'(x) 

= E'^(x) 


( 6 ) 


■ Ai' 
B^ 


fj(E,A,</,o) 


r 4™ 1 


(15) 


where v(x) = vf/(1 + ^i(x)) is the renormalized Fermi 
velocity. For constant A, there are two independent 
eigenstates with energy E given by 


q,±(x) V^{E,A) 


w- 


'(x) = f du 

Jo 


V{u-ao)±y/E‘^-\A\^ (7) 

v{u) 


77 +(A, A) 


I 


77 -(A, A) 


I 


yjE + y/W^W 
-A/yjE+y/W^W _ 

-A* /yjE+yEE^-\A \^ 

yjE+,/W^W 


To model the Josephson Junction, we consider the 
stepwise variation of A(a:) in three regions a: < 0 (re¬ 
gion I), 0 < a: < ccp (region II) and x > Xp (region III) 
given by 


A(a:) = 


( 10 ) 


where the wave function on the two ends of the 
Josephson junctions are connected by Tj{E, A, (j)Q) = 
TBiE, Ao)-^ fnixp)-^ Tb{E, Aoe*’^'>). The S'-matrix for 
the junction must relate the incoming and outgoing states 


as 


Aiii 

B^ 


= S 


Ai 

5111 


; this 5'-matrix can be written in 


terms of Tj{E, A, ^ 0 ) as 


5 = 


Tj{E, A, (()o)ii 


1 

Tj{E, A, 4 ) 0)21 


—Tj{E, A, 4>o)i2 


(16) 

The supercurrent in the junction can be derived us¬ 
ing the well-known relationship between the Josephson 
current of the junction and the spectrum [^, namely. 


I = h+l 2 +l 3 (17) 

(EpI^kBT) ^ (18) 

h = -l r dE ln[2cosh(F;/2fcBT)]^^^l^(I9) 
^ Jao d<l> 

h = J dx\A{x)\^/g 


Ao (1) 

0 (II) 

Aoe^^o (III) 


( 20 ) 
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/i is the contribution form the discrete spectrum of in¬ 
gap states and I 2 is the contribution form continuum of 
states with energy above the gap with density of states 
p{E, (j)) for the one spin state at each fermi point as we 
have for the helical metal. The third them I 3 vanishes 
for the phase independent gap and will be ignored in 
this letter. In (l20)l . g is the interaction constant of BCS 
theory. 

For the states with energy E < Aq, k is imaginary and 
as a result 'I'^(a;) in region I and in III are not 

normalizable. Equation (|15p then simplifies to 


0 


= TsiE, Ao)-^fnixp)-^ TsiE, Aoe'^'>) 


^iii 

0 

( 21 ) 

and leads to the following equation determining the in¬ 
gap energies: 


En 


(t>0 


cos ^(F;„/Ao)-I -^£-|-^=n7r, 

V 2 


n G Z (22) 


where C = fj’’ dx (1 -I- Ei {x)) is the effective length of 
the junction which is modified by fluctuations of the im¬ 
purity. This is the new “optical length" of the junction. 
The phase factor (j)i has cancelled out in confirm¬ 
ing that static impurities have no effect on the energy of 
in-gap Andreev states. The effect of impurities is only 
through their dynamical fluctuations that leads to the fi¬ 
nite self-energy Ei(a:) which modifies the effective Fermi 
velocity. The effect on the energy eigenvalues is most 
vividly illustrated by considering states with ^ Aq, 
in which case we get En ~ {v/2C) [27r(n + ^) — (Ao] • The 
increase in C implies that En and dEn/dcj^o are decreased 
relative to the case with no impurities. More generally, 
defining 0„ = + cj)oj 2 — rnr, the supercurrent asso¬ 

ciated with each in-gap state reads as 


Tn(0o) = - 


edEn e|Ao| 


h 


2h 


sin(0„) 


1 + sin(0„)£|Ao|/?; 


As a function of 0^, , this has a maximum at 0, 


(23) 

= 7r/2, 

so that the critical current is Icrzi = (e|Ao|/S)(l -I- 
AlAol/f)”^. The condition 0„ Ri 7r/2 is actually ob¬ 
tained for modes of very low energy En <C |Ao|. It is 
important to note that by that the suppercurrent gen¬ 
erated by in gap states decreases by increasing C which 
shows that impurities clearly affect the supercurrent. 

For the states above the gap, the density of states 
is given by the Krein-Friedel-Lloyd formula p{E) = 
(lndet5) (s^. Using (fTBll we get 

1 - /3|cos^(0£;) 


det5 = e-*‘^“- 


j3% cos(0b) - zsin(0_E) 


where 0 _e = ^A-|-0o/2 and Pe = ^- The supercurrent 
due to above-the-gap states then simplifies to 

^ \r iE,.nHEI2k,T)2A^ 


E = 


eksT 

2Trh 


In [2 cosh (Ao/2fcBT)] 


e r°° 

-r^ dE tanh{E/2kBT) 
dTTft 


I Aq 

eksT 

2TTh 


E^E^ - Ag 
E "^ — Aq cos^(0b) 


In [2 cosh (Ao/2fcBr)] 


- 1 


(24) 


We would like to emphasize two important features 
of the supercurrent contribution from states with energy 
above the superconducting gap: 

1. For low temperatures T ^ Aq, I 2 is only weakly 
T-dependent through the temperature dependence 
of superconducting gap Aq. 

2. I 2 is also only weakly dependent on A, i.e, only 
weakly sensitive to impurities. 

To elucidate the second point, we first note that the 
second term in (IMl) is independent of C. Assuming 
T <g; Aq, tanh (i5/2fcBT) ~ 1. The integrand in the 
first term in (l24)) has two types of dependence on E. 
One is a periodic dependence, with period hv/C due 
to cos(0_e), and the other is a decaying dependence, of 
the form for large E. For the effective junction 

length C larger than Hv/Aq, the oscillatory dependence 
is is much faster than the decay rate and so can be aver¬ 
aged over 0 E ■ (This may be viewed as an application of 
the Riemann-Lebesgue lemma.) As a result, the depen¬ 
dence on C will be eliminated and I 2 will not be seriously 
affected by impurities even when fluctuation effects are 
included. 

In conclusion we have shown that the supercurrent in 
Josephson junctions with helical metals, such as on the 
surface of three-dimensional TIs, is affected by impuri¬ 
ties through their temporal fluctuations. However, this 
applies primarily to the supercurrent generated by in¬ 
gap Andreev states. The supercurrent carried by the 
states above the gap will not be seriously affected by 
impurities. Based on our results, the supercurrent in 
the Josephson junctions on the surface of TIs can be in¬ 
terpreted as a superposition of two contributions, one 
which is strongly temperature-dependent and also sen¬ 
sitive to the impurities in the junction and one which 
is only weakly temperature-dependent and not sensitive 
to the impurities. Given new advances in controlling the 
level of disorder in TIs [Hill , these results will be useful 
in analyzing many of the experimental results on Joseph¬ 
son junctions made on TIs. For example, our analysis 
is consistent with the experimental results in [2lJ, l22l |: 
whether different level of impurities could affect the crit¬ 
ical current in the Josephson junction on TI was the main 
missing ingredient in the theoretical model used to inter¬ 
pret those results. In fact, our work may be considered 
as further substantiating the interpretation, presented in 
2l| . in terms of two types of supercurrent contributions. 
The detailed comparison with those experimental results 
will be subject of a following publication. 
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Supplementary Materials 


I. EFFECTIVE ACTION FOR ELECTRONS FROM OSCILLATING ELECTROSTATIC IMPURITY 

A. The perturbative result 


Here we present the calculation of the effective action for the electron field due to the interaction with impurities. 
The starting action is the one given in Eq.(3) in text. Since we will also be discussing how temperature affects the 
calculation, it is more convenient to use a Euclidean space (imaginary time) formalism, so that thermal effects can be 
included via the Matsubara formalism. 

The Euclidean version of (3) is given by 


S' = y cPx (T>4 — i T3D1) ^ + J d 

Defining 7 ° = 7 ^ = ti, 7 ^ = —T 2 , this becomes 


2 ™ 


0 A 
A* 0 


<i> 


S = Jd^x^ [7^(9^ - fra A^) + M] $ 


where $ = $ 17 ° and 


M = 


A* 0 
0 A 


The propagator for the electron is then given by 
S{x,y) = {<P{x)^{y)) = 


j-d + M 


f _(PP_ (f7 ■ P + M*) 


(2r)2 


(p 2 + M*M) 


(SI) 


(S2) 


(S3) 


(S4) 


The electrostatic field can be expanded as ^ 4 ( 0 ; + ^) ~ A 4 {x) + ^ 51^4 = A 4 {x) + ^ F 14 . The correction to the action, 
to quadratic order is then given by the Wick contractions of —^SintSint', we need one electron propagator and one 
contraction for the f’s. Thus 


^Seff = 


Fi4(a;)Fi4(y) (^(a;)^(?/)) d>(a;)7°T3 S(a;, ?/)7°T3d>(y) 


(S5) 


For us, ^ depends only on time as it is the oscillating coordinate of the impurity relative to the mean position. Further, 
if we consider several impurity atoms, only the ^’s of the same impurity atom can have nonzero average {^ix)^{y)) . 
This means that the contribution to the integral is concentrated around Xi = yi- We van encode these by wriitng 


(C(a;)^(y)) = 


E 


/ 


drk 


^-ik-{x-y) _ 


1 


Mi J (27r)2 


fc? + ( 


(S 6 ) 


The spatial part gives a delta function and integrating over xi, we get the factor R. This shows that R may be 
taken as a rough measure of the extent over which this correlation exists. Using this in (|S5ll and changing variables 
p = q — k, we get 

^Seff = I 7 ^e-*« ("-^)$(i)E(g)$(y) (S7) 


S(g) = 


( 2^)2 

e^RF^^ f d?k 7 ^X 3 (*7 • (g - fc) + M*) 7 ^X 3 
Mi J (27r)2 {kj + uj^)[{q - k)^ + M*M] 
e^RF ^4 


f 

Jo 


da 


d^K (f 7 °g 4 — M){1 — a) — aM 


(5 8 ) 

(59) 


Mi Jo ““ J (27r)2 [A| + aKf + Q2]2 

= a(l — a)ql + (1 — a)uj‘^ + aM* M 

In the last line of (ISSp we have combined the denominators of the previous expression using Feynman’s formula and 
also defined K 4 = k 4 — aq 4 , Ki = ki — qi. Notice that the self-energy S is naturally split as 

E = -El (f7°g4 - M) -k E 2 M 
















2 


51 

5 2 


e'^RF^^ f d^K (1 - a) 

Mi Jo J (271)2 [^2 ^ + q 2]2 

e'^RF^i r d'^K a 

Mi Jo J (271)2 [^2 ^ ^^2 + q 2]2 


(510) 

(511) 


Evidently, the small q limit of Ei will give a term of the form $( 7^94 + M)$ in the action, while the same limit of E 2 
will correct the mass M. Higher order powers of qi in the expansion of E will give higher derivatives of the electron 
held; these are not important since we are interested in the low energy modes of the electron. Explicitly, 


ASeff ^ Ei(g 4 = 0 ) 


$(7°a4 + M)$ +E2(g4 = 0) 




(S12) 


Carrying out the K 4 integration. 


Ei(g4 = 0) 



e^F^pR 1 

Mi 27iw2 


a) 


dKi 1 

871 [aKf + aM*M + (1 — Q:)a;2]i 


(S13) 


In the last line, we evaluated the remaining integral neglecting M*M in comparison to uj'^, since the frequency of 
vibration for impurity atom is much larger than the the possible gap. In any case, we set M =) for the region II of our 
discussion in text to which this calculation is applied. We also made the continuation to real time by F^^ —>■ F^p = 

The result (IS13I) agrees with what is quoted in text. 

In a similar way. 




(S14) 


We can include hnite temperature effects by using Matsubara frequencies for = 27rri//3 = a;„ and p 4 = 27r(n + 
P = I/(fcsT). The result is then 


^ E / E(g) $(y) 

_ e^Ff^R 1 r dki 77°(g4^ - a;„) + M* + 

^ Mi P ^ J 27r (u;2+a;2)[(g4s - w„)2 + + M*M] 


(S15) 


As usual, we carry out the summation via contour integration. Writing = kf + M*M, 

1 *7°(94s - w„) + M* + 77^fci _ 1 ^ 77°(g4s - z) + + M* 

P ^ {ujI+ w2)[(g4s - a;„)2 + kl + M*M] 27r / - I) (2;2 + a;2)[(g4s - z)2 + Q 2 ] 


where the contour encloses the poles of the (e*^^ — 1)“^ factor. Folding the contours back to the upper and lower 
half-planes and evaluating the residues at the other poles, we find some terms which are independent of P (and coincide 
with what we have already done in (IS 8 I IS14|) 1 and a set of terms which are /^-dependent. The /3-dependent part of 
El is 


1 _ e^Ff^R r dki 

]\/[j j 

The key point is that, for us, lu ^ ksT. As a result, the T-dependent correction is exponentially suppressed due 
to the factor in the denominator. A similar argument holds for E 2 as well. Thus, in conclusion, the effective 
action is of the form as in Eq.(5) of the text where we can take Ei and E 2 to be independent of temperature. Any 
temperature-dependence of the supercurrent would be due to the T-dependence of the gap A and due to factors such 
as tanh(i3/2fcBr). 


+ (<74s - *w)2 <52 _|_ _|_ j^)2 


o/Sw _ 


(S17) 
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B. The more general argument 


The form of the action in Eq.(5), which is all we need for the rest of the results in the paper, can be obtained 
on general symmetry grounds. We will first consider a region of uniform distribution of impurities. The starting 
action has a (1 + l)-dimensional Lorentz symmetry with the Fermi velocity in place of the speed of light. If the 
interactions respect this symmetry, the relative coefficients of the two terms in the combination + ^T 2 dx^ 

must be preserved. However, the interactions do not preserve this Lorentz symmetry and so, on general symmetry 
grounds, we expect the effective action to be of the form 


Seff = A / + B / 


/■ 




+terms of higher order in derivatives of $ 


(S18) 


The higher derivative terms are not important for the low energy modes which are of interest to us. H, H, C are 
calculable constants. We can scale out one of them, say, B to write 


— 





+ 




(S19) 


We now define Ei 


H/H = l + Ei 


(S20) 


With this, we see that the first term (IS19I) has the form in Eq.(5). 

As for the term, notice that the derivative terms of the starting action, Eq. (3) in text, has a chiral 

symmetry, $ —>■ This is true even with the electromagnetic interactions. Thus, if M is originally zero, it 

cannot be generated by perturbative corrections. Therefore, in the theory with nonzero, C must be such that it 
vanishes when M —>■ 0. We therefore write C ~ M, and define 112 in general by 


C 

B 


(1 + Si)M + MS2 


With the two definitions (IS20II and (jS21(l . we get 


^eff — 


j dxdt [1 + Tji{x)) ^ (^Tidt + M {x)^ 


$ 


+d>T2i9a;<i) + $ M(a;)E2(a;) d) 


(S21) 


(S22) 


in agreement with what is given in text. The argument of the previous subsection was given to see how the perturbative 
calculation of the coefficients A/B and C/B can be done and give results consistent with the general expectations. 










